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\ An analogy is developed between the breaking of space-inversion symmetry in equilibrium statis- 

tical mechanics and the breaking of time-reversal symmetry in nonequilibrium statistical mechanics. 
In this way, similar relationships characterizing fluctuations are obtained in both contexts. 
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I. INTRODUCTION 



Symmetry breaking is a major phenomenon in natural systems. A priori, the equations describing a system may 
be symmetric under a group of continuous or discrete transformations. For some reasons, the system is often found in 
some state that is a solution of the describing equations without satisfying all the symmetries of these equations. In 
such cases, the symmetry is broken by the particular state of the system. Symmetry breaking phenomena may have 
different origins such as external forces or constraints imposed to the bulk or the boundaries of the system. In other 
cases, the symmetry may be broken by the initial conditions that are selected. In spontaneous symmetry breaking, 
the symmetric state is unstable so that the system has to evolve towards some stable but asymmetric state, which is 
\ selected by slightly biased external perturbations or initial conditions. 

. These different phenomena can be studied for the symmetries that are broken. They can be continuous leading to 
, ' ' the emergence of Goldstone-Nambu massless modes (e.g., the transverse acoustic modes in solids) or Anderson-Brout- 
Englert-Higgs massive modes if a local gauge symmetry is broken in the presence of the corresponding long-range 
£h ■ interaction. 

Here, we focus on discrete Z2 symmetries such that the symmetry group is composed of the identity 1 and an 
involution X. Z 2 = {1, X}. Our purpose is to compare the breaking of time- reversal symmetry by thermodynamic 
forces in nonequilibrium steady states to the breaking of space inversion by external fields in equilibrium states. In 
. the later case, we aim at obtaining fluctuation relations similar to those known for nonequilibrium steady states, using 
the analogy between them through the phenomenon of symmetry breaking. 

The paper is organized as follows. In Section [Til the analogy between broken space-inversion and time-reversal 
symmetries is illustrated in simple Hamiltonian systems with one degree of freedom or infinitely many noninteracting 
particles. In Section IlIIl fluctuation relations are deduced for equilibrium spin systems in an external magnetic field 
breaking a space-inversion symmetry. In Section IIV1 nonequilibrium fluctuation relations are presented in analogy 
l/^ I with the previous ones, emphasizing the role of the time-reversal symmetry breaking by the invariant probability 
measure describing the nonequilibrium steady state. Conclusions and perspectives are drawn in Section IVl 

O 

(N) ' II. SYMMETRY BREAKING BY SELECTION OF INITIAL CONDITIONS 

In this section, we start by considering two Hamiltonian systems with one degree of freedom illustrating Z2 symmetry 
breaking. The first is the Hamiltonian flow in a double-well potential and the second is the simple pendulum. Both 
Hamiltonian systems are symmetric under space inversion x — > —x (i.e., parity) and time reversal t — > —t. However, 
there exist orbits breaking the symmetry under parity in the first system and the symmetry under time reversal in 
the second system. As a consequence, the invariant probability measures supported by these orbits also break the 
corresponding symmetry. In these one-dimensional systems, the orbits are periodic except the separatrices formed by 
the homoclinic orbits connected to the unstable stationary points. For the periodic orbits, an invariant measure can 
be defined by temporal averaging over the period T to get the invariant probability density 

p c (r) = ±£ 5(r-&r )dt (i) 

where T — (x,p) denotes a point in phase space, Tq is the initial condition of the orbit, and $' is the Hamiltonian 
flow. The support of the so-defined probability measure is the curve depicting the periodic orbit in the phase space: 



c = {r = $ t r : te [o,r[} 



(2) 
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We may wonder if the invariant density ([T]) coincides with the microcanonical invariant density defined by 

S [E - H(T)\ 



Pe(T) 



jS[E-H(T)]dT 



(3) 



where H(T) = H(x,p) is the Hamiltonian function and E = H(Tq) is the energy fixed by the initial condition Tq. 
Let us look at this question in the two aforementioned examples. 

1. Double- well potential. For this system, the Hamiltonian function is given by 

a n b A 
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U(x) 



where 



U(x) 



■ x 



(4) 



with a > and b > 0. The phase space is the plane M. — {T = (x,p) £ M 2 }. The Hamiltonian function is invariant 
under the parity transformation: 



n(a;,p) = (-x, -p) 



(5) 



which is an involution since n 2 = 1, generating the group Z2 = {l,n} of symmetries. Consequently, Hamilton's 
equations of motion stay invariant under parity. 

The phase portrait of this system is plotted in Fig.QJ,. The system has three stationary points: the origin x = p = 
which is symmetric but unstable and the two minima of the potential at x = ±yj a/b and p — 0, which are stable 
at the minimum energy U m - ln = —a 2 /(4b). At zero energy, there are two homoclinic loops connected to the unstable 
point and they play the role of separatrices between the periodic orbits at negative and positive energies. 

With every positive value of the energy E > 0, there corresponds one and only one periodic orbit Ce that oscillates 
between both wells. In contrast, two periodic orbits correspond with every negative value of the energy C/ m i n < 
E < and they oscillate in either the well on the right-hand side or the one on the left-hand side. If the periodic 
orbits at positive energies are symmetric under parity, this is no longer the case for those at negative energies: 



ILC E = C E 
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(G) 
(7) 



Therefore, every periodic orbit at negative energy breaks the symmetry under space inversion of Hamilton's equations. 





(a) (b) 
FIG. 1: Phase portraits of the Hamiltonian flows for (a) the double- well potential ((4]) and (b) the simple pendulum 



2. Simple pendulum. The Hamiltonian function is here given by 

P 2 

H = — mgl cos x (8) 

2ml 2 

where m is the mass of the body attached to a massless rod of length / in the gravitational acceleration g oriented 
towards x = 0. The angle varies in the interval — ir < x < +ir so that the phase space is the cylinder M. = {T = (x,p) : 
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x g] — 7r, +tt],p G R}. The Hamiltonian function and Hamilton's equations are symmetric under the time-reversal 
transformation: 

Q(x,p) = (x,-p) (9) 

which is also an involution since O 2 = 1, generating the group Z2 = {1, 0} of symmetries. 

The phase portrait of the pendulum is plotted in Fig. [T]d. The system has two stationary configurations: the 
stable one x = p = at the minimum energy U m - m = —mgl and the unstable one (x = ±7r,p = 0) at the energy 
tAnax = +mgl. On the cylinder, the two unstable points {x = ±7r,p = 0) are equivalent so that the two curves 
connecting them form two homoclinic orbits. They constitute the separatrices between the oscillating periodic orbits 
at lower energies U m i n < E < C/ max and the rotating periodic orbits at higher energies E > C/ max . 

With every value of the energy i7 m i n < E < C/ m ax, there corresponds one and only one oscillating periodic orbit 
Ce which is self-reversed because the time-reversal transformation maps it onto itself. In contrast, there exist two 
periodic orbits rotating either anticlockwise or clockwise corresponding with every energy value E > C/ max . None 
of them is self-reversed so that each one breaks the time-reversal symmetry of Hamilton's equation. In analogy with 
Eqs. ©-0, we have that 

QC E = C E for C/ min < E < U max , (10) 

9C4 +) =C E - ] ^C { E +) for E>U max . (11) 

The simplicity of this example shows that the breaking of time-reversal symmetry is a common phenomenon arising 
by the selection of initial conditions. The fact is that the symmetry of the equations of motion does not imply the 
symmetry of all their solutions. 

As a consequence, the invariant probability density ([T]) also breaks the symmetry if it is supported by one of the 
periodic orbits C E : 

p&e+) (er) = p C (_) (T) + p c (+, (T) (12) 

for E > U max . More generally, it is possible to define a continuous family of invariant densities supported by both 
periodic orbits according to 

p A = A p c (+) + (1 — A) p c (- } (13) 
with < A < 1. They satisfy the symmetry relation 

PA (er) =pi_ A (r) (w) 

and they break the time-reversal symmetry unless A = 1/2. In this latter case, we recover the microcanonical invariant 
density (J3J), which is always symmetric under time reversal since it is defined in terms of the symmetric Hamiltonian 
function: px=i/2 = Pe = Pe ° 6. 

3. Effusion. In order to show that the previous considerations extend to situations relevant to nonequilibrium 
statistical mechanics, let us envisage the effusion process of an ideal gas of noninteracting particles through a small 
hole in a thin wall separating two infinite domains in the three-dimensional physical space [l[ . The process is depicted 
in Fig.fSJ In the domain on the left-hand side, the particles arrive from infinity with velocities and positions distributed 
according to a Maxwell-Boltzmann distribution at the temperature Tl and particle density n^. Most of these particles 
undergo specular collisions on the wall and return to infinity in the same domain. However, some particles cross the wall 
through the hole towards the other domain. Vice versa, the domain on the right-hand side contains particles coming 
from the opposite infinity with velocities and positions distributed at the temperature Tr and particle density tir. 
Most of them remain in the same domain after one specular collision on the wall, while some of them cross the wall 
through the hole. 

The motion of every particle is ruled by the Hamiltonian flow $* in the single-particle phase space M. = {T = 
(r,p) G (R 3 — W) ® M 3 } where W denotes the thin wall with its hole. The Hamiltonian function is given by 
H = p 2 / (2m), but the orbits undergo specular collisions on the wall. Accordingly, the Hamiltonian flow is symplectic 
and time-reversal symmetric: 8 o f ' o 8 = $ _t . 

The single-particle distribution function has the Maxwell-Boltzmann form: 

n c _ ( P 2 \ 



« r >=^>=<w^ e n-3^J (15) 
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FIG. 2: Schematical representation of the effusion process in the physical position space r = (x,y,z) £ R . The wall W 
separates the position space in two domains where the gas of noninteracting particles has different temperatures and densities. 
The figure depicts three types of phase-space orbits d and their time reversal Qd (with i — 1 , 2,3). They are here projected 
onto the position space by the projection P : (r, p) € R 6 -> r 6 R 3 . 



where the particle density nc and the temperature Tc take the values associated with the orbit C to which the phase- 
space point r = (r, p) belongs and corresponding to the domain from which the orbit is coming, (fee is Boltzmann's 
constant.) For the three types of possible orbits depicted in Fig. ® these values are the following: 

n C = n L , T c = T L for r e d , ®C X , C 2 ; (16) 
n c = n R , T c = T R for T E C 3 , 0C 3 , QC 2 . (17) 

The distribution function is invariant under the Hamiltonian flow, /($*T) = f(T). Nevertheless, it is not symmetric 
under time reversal, /(Or) ^ /(r), for nonequilibrium constraints such that ^ 7Jr or Tl ^ Tr. 

For the infinite-particle system, an invariant probability measure /z can be constructed as a Poisson suspension on 
the basis of the invariant distribution function (|15p , as explained in Ref. [2j . The so-constructed infinite dynamical 
system ($^ , A4°° , fi) is not symmetric under time reversal because /jo0 / (i under nonequilibrium constraints, 
although the flow itself always is, as pointed out in Ref. [||. 



III. SYMMETRY RELATIONS FOR EQUILIBRIUM STEADY STATES 

In this section, we consider spin systems at equilibrium as described by the Hamiltonian function: 

H B (a) = H (a) - BM N (a) (18) 

where a = (ui,(J2, ■■■,&n) is a configuration of the spins <jj — ±1, Ho(a) is a Hamiltonian that is symmetric under 
the group Z2 = {l,n} with Tla — —a and another group G of spin permutations including translations of the type 
T{<J\, 02, C3..., (Ta?) = (ctat, (Ti, (72-, <7/v-i)- Moreover, B denotes an external and uniform magnetic field and 

N 

M N (<r) = J2<T3 (19) 
3=1 

is the total magnetization. Accordingly, the total Hamiltonian Hb(ct) is symmetric under the group G, but the Z2 
symmetry is broken if B ^ 0. We suppose that the system is in the Gibbsian equilibrium state 

MW = ^e^W (20 ) 

with the inverse temperature /3 = (/cbT) -1 and the partition function Z — tr q~0 h b{ o ) _ This equilibrium state has 
the same symmetries as the Hamiltonian. The fluctuations of the magnetization ((19)) can be characterized in terms 
of large-deviation functions [3, Q . 
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The main result is that, in the thermodynamic limit N — > oo, the equilibrium probability distribution Pb(M) of 
the magnetization (|19p obeys the symmetry relation: 

MM) ^ 20BM . 



Pb(-M) 

Introducing the generating function of the magnetization cumulants as 

Q B {\)= lim -lln(e- AM «) B (22) 
N— >oo Jy 

the symmetry relation (1211) yields the equivalent form: 

Qb(X) = Q B {2pB-X) (23) 

because the large-deviation functions Pb(M) and Qb(A) are the Legendre transforms of each other. Equations (|2"Tj) 
and (l23l) express the breaking of the Z2 symmetry by the external magnetic field. These results are established here 
below for two different models by using their large-deviation properties reported in Ref. [J] . 

1. Ising model. The spins form a one-dimensional chain and their interaction energy is given by 

N N 

H B (a) = -jJ2 *J o-i+i -BY,** ( 24 ) 
3=1 3=1 

with CT/v+i = (Ti. This Hamiltonian function is symmetric under the group of translations G = {1, T, T 2 , r N_1 }. 
Now, the generating function (|22[) can be calculated in terms of the largest eigenvalue of the transfer matrix 

/ Q pj+f3B-X e -/3J \ 
V * = ( e -PJ e /3J-/3B+\ ) ( 25 ) 



which obeys the symmetry relation 



7r V\ 7T = V2/3B-A with i = ( J J ) . (26) 



Accordingly, its eigenvalues have the symmetry A — > 2/3B — A and the relation ([23]) is satisfied. 
2. Curie- Weiss model. Here, the Hamiltonian function can be written as 

H B (tr) = -^M N (a) 2 - B M N (a) (27) 

in terms of the magnetization (I19[) . This Hamiltonian is invariant under the symmetric group G — SymTV, which 
contains subgroups similar to the translation group. Moreover, the Z2 symmetry is broken if B ^ 0. The number of 
spin configurations such that the magnetization takes the value Mm = N — 2n is equal to g n = N\/[n\(N — n)\] with 
n = 0, 1, 2, N. In the thermodynamic limit N —> 00, the generating function (|22[) is given by 

Qb(A)= hm lln gjg, (28) 

with 



,+1 

W^jvCO = J dm exp 



, . . P J 2 1 ~ m 1 — 771 1 + m 1 

iv 777 + t 777 In In — 

"2 s 2 2 2 



(29) 



where 777 = 1 — 2n/N stands for the magnetization per spin. Now, the symmetry Wn(£) = Wn{— £) of the partition 
function (|29[) implies the symmetry relation (|23j) . hence Eq. (|21[) . Spontaneous symmetry breaking only happens 
if _B = and below the critical temperature of the paramagnetic-ferromagnetic transition where the magnetization 
probability distribution Pg(M) is bimodal. 
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IV. SYMMETRY RELATIONS FOR NONEQUILIBRIUM STEADY STATES 

The symmetry relation (|2ip is analogue to relationships established for nonequilibrium steady states such as the 
current fluctuation relation: 

^ oc exp(A.Ji) (30) 

which has been proved for all the fluctuating currents J flowing across some mesoscopic device under nonequilibrium 
constraints fixed by the thermodynamic forces, also called the affinities A [(J. In terms of the cumulant generating 
function 

Qa(A) = lim -iln(c- A J, )A (31) 

t— too t 

the relation (1301) implies the equivalent expression: 

Qa(A) = Q a (A-A) (32) 

which is similar to Eq. (f2"3"f . Here, the relations (I3U1) and (|3"2")l express the breaking of the time-reversal symmetry in 
the nonequilibrium steady state driven by non- vanishing affinities A. The principle of detailed balancing is recovered 
at equilibrium for A = 0. These relations have fundamental consequences on the linear and nonlinear response 
properties Q. An example is given by the effusion process with the invariant probability measure constructed here 
above, as shown in Ref. \M- 



V. CONCLUSIONS AND PERSPECTIVES 



In this paper, the analogy has been developed between different kinds of broken Z2 symmetries. In Section [TTl 
simple Hamiltonian flows already illustrate these symmetry-breaking phenomena induced by the selection of initial 
conditions and their profound similarities. In Sections Mil and IIV[ this analogy has been further emphasized. On the 
one hand, Eqs. (1211) and (|23l) have been established for two equilibrium models where the space-inversion symmetry 
is broken by an external magnetic field. On the other hand, the similar relationships ([30)) and (|32|) characterize the 
current fluctuations in nonequilibrium steady states where the time-reversal symmetry is broken by some external 
drivings A. These results express the effects of the broken Z2 symmetry on the fluctuations of the relevant quantities, 
showing the generality of such considerations, which concern other Z2 symmetries as well [8]. 
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